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Recent results from the KASCADE experiment on measurements of cosmic rays in the energy range of the 
knee are presented. Emphasis is placed on energy spectra of individual mass groups as obtained from an two- 
dimensional unfolding applied to the reconstructed electron and truncated muon numbers of each individual 
EAS. The data show a knee-like structure in the energy spectra of light primaries (p, He, C) and an increasing 
dominance of heavy ones (A <; 20) towards higher energies. This basic result is robust against uncertainties of the 
applied interaction models QGSJET and SIBYLL which are used in the shower simulations to analyse the data. 
Slight differences observed between experimental data and EAS simulations provide important clues for further 
improvements of the interaction models. The data are complemented by new limits on global anisotropies in 
the arrival directions of CRs and by upper limits on point sources. Astrophysical implications for discriminating 
models of maximum acceleration energy vs galactic diffusion/drift models of the knee are discussed based on this 
data. 



1. Introduction 

A puzzling and most prominent feature of the 
cosmic ray (CR) spectrum is the so-called knee, 
where the spectral index of the all-particle power- 
law spectrum changes from approximately —2.7 
to —3.1. Several models have been proposed in 
order to explain this feature shown in Fig.^ but 
none of them has managed to become broadly ac- 
cepted. Some models focus on a possible change 
in the acceleration mechanism at the knee |1I2I3| , 
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e.g. due to the limiting energy defined by the size 
and magnetic field strength of the acceleration re- 
gion (E max Si Zx(BxL)). Other ones discuss an 
increased leakage of CRs from the Galaxy due to 
a change in the confinement efficiency by galac- 
tic magnetic fields, e.g. |4I5I6I7| . Again, this re- 
sults in a rigidity scaling of the knee according 
to the maximum confinement energy. Finally, a 
third group of models attributes the effect of the 
knee to CR interactions at their sources, during 
their propagation in the Galaxy, or in the up- 
per atmosphere. Such scenarios include nuclear 
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Figure 1. Compilation of the all-particle cosmic 
ray spectrum showing the knee, the suggested 
second knee, and the ankle of the CR spectrum 
(compilation by H. Ulrich). 



photo-disintegration processes by UV-photons at 
the sources jS], interactions of CRs in dense fields 
of massive relic neutrinos [9], production of gravi- 
tons in high-energy pp collisions etc. A re- 
cent review about this topic can be found e.g. in 
Refs. P3E21- 

To distinguish between these models and allow- 
ing to answer the long pressing question about the 
origin of cosmic rays and about the knee in their 
spectrum, high quality and high statistics data 
are required over an energy interval ranging from 
at least 0.5 to 500 PeV. It appears worthwhile to 
mention that solving the old problem about the 
origin of CRs in the PeV region is a prerequi- 
site also for an understanding of the highest en- 
ergies in the GZK-region. Due to the low flux 
involved at energies <; 10 14 eV, only extensive 
air shower (EAS) experiments are able to provide 
such data. In EAS experiments, primary CRs are 
only indirectly observed via their secondaries gen- 
erated in the atmosphere. The most important 
experimental observables at ground are then the 
electromagnetic (electrons and photons), muonic, 
and hadronic components. In addition or alterna- 
tively, some experiments also detect photons orig- 
inating from Cherenkov and/or fluorescence radi- 
ation of charged particles in the atmosphere. For 
a brief review about EAS observables and their 
experimental techniques the reader is referred to 



Refs. puma. 

Unfortunately, progress on interpreting EAS 
data has been modest mostly because of two rea- 
sons: Firstly, the EAS development is driven both 
by the poorly known high-energy hadronic inter- 
actions and their particle production in the very 
forward kinematical region as well as by uncer- 
tainties in the low energy interaction models influ- 
encing mostly the lateral particle density distribu- 
tion functions (TJ|. Secondly, due to the stochas- 
tic nature of particle interactions, most impor- 
tantly the height of the very first interaction in 
the atmosphere, EAS are subject to large fluc- 
tuations in particle numbers at ground. To make 
things even more complicated, the amount of fluc- 
tuations depends, amongst others, sensitively on 
the primary CR energy and mass jl^j . Here, it is 
very important to realize that EAS fluctuations 
are not to be mistaken as random Gaussian er- 
rors associated with the statistics in the number 
of particles observed at ground. The latter one 
can be improved by the sampling area of an EAS 
experiment, while the former one is intrinsic to 
the EAS itself, carrying - for a sample of events 
- important information about the nature of the 
primary particle. Clearly, both kinds of fluctua- 
tions have to be accounted for in the data analysis 
of steeply falling energy spectra in order to avoid 
misinterpretations of the observations. 

2. Results from the KASCADE Experi- 
ment 

KASCADE (Karlsruhe Shower Core and Array 
Detector) is a sophisticated EAS experiment for 
detailed investigations of primary CRs in the en- 
ergy range of the knee. For reconstructing the 
CR energy and mass and for investigating high- 
energy hadronic interactions, KASCADE follows 
the concept of a multi-detector set-up providing 
as much complementary information as possible 
as well as redundancy for consistency tests. Most 
relevant for the results presented in this paper is 
the scintillator array comprising 252 detector sta- 
tions of electron and muon counters arranged on 
a grid of 200 x 200 m 2 . In total, it provides about 
500 m 2 of e/7- and 620 m 2 of //-detector coverage. 
The detection thresholds for vertical incidence are 
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Ef* ~ 5 MeV and Ef T ~ 230 MeV. More details 
about the e/7- and /x-detectors and all other de- 
tector components can be found in Ref. |1K| . 

2.1. Chemical Composition and Energy 
Spectra 

The traditional and perhaps most sensitive 
technique to infer the CR composition from EAS 
data is based on measurements of the electron 
(N e ) and muon numbers (N^) at ground. It is 
well known |13j that for given energy, primary Fe- 
nuclei result in more muons and fewer electrons at 
ground as compared to proton primaries. Specifi- 
cally, in the energy range and at the atmospheric 
depth of KASCADE, a Fe-primary yields about 
30 % more muons and almost a factor of two 
fewer electrons as compared to a proton primary. 
The basic quantitative procedure of KASCADE 
for obtaining the energy and mass of the cosmic 
rays is a technique of unfolding the observed two- 
dimensional electron vs truncated muon number 
spectrum of Fig.|2]into the energy spectra of pri- 
mary mass groups [17) . The problem can be con- 
sidered a system of coupled Fredholm integral 
equations of the form 

+00 

dJ f (Ma 

d\gN e dlgiVtr =2^ J d l gE 

A -00 

• PA QgNe, lgiV* r I lgE)-d\gE 
where the probability pa 
PA (]gN e ,]gN^ I IgE) = 

—00 

is another integral equation with the kernel func- 
tion kA = ta ■ (-a • S A factorizing into three parts. 
Here, ta describes the shower fluctuations, i.e. 
the 2-dim distribution of electron and truncated 
muon number for fixed primary energy and mass, 
tA describes the trigger efficiency of the exper- 
iment, and sa the reconstruction probabilities, 
i.e. the distribution of N e and A^. t , r that is recon- 
structed for given true numbers JV|, Nff' 1 of elec- 
tron and truncated muon numbers. The probabil- 
ities pa are obtained by parameterizations of EAS 
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Figure 2. Two dimensional electron (N e ) and 
truncated muon number (iV* r = J* 40m Pix{ r )dr) 
spectrum measured by the KASCADE array. 
Lines display the most probable values ex- 
pected for proton and iron primaries according 
to CORSIKA simulations employing two differ- 
ent hadronic interaction models [17) . 



Monte Carlo simulations for fixed energies using 
a moderate thinning procedure as well as smaller 
samples of fully simulated showers for the input of 
the detector simulations. Because of the shower 
fluctuations mentioned above, unfolding of all 26 
energy spectra ranging from protons to Fc-nuclci 
is clearly impossible. Therefore, 5 elements (p, 
He, C, Si, Fe) were chosen as representatives for 
the entire distribution. More mass groups do not 
improve the ^-uncertainties of the unfolding but 
may result in mutual systematic biases of the re- 
constructed spectra [P7] . 

The unfolding procedure is tested by using ran- 
dom initial spectra generated by Monte Carlo 
simulations. It has been shown ^7] that knee 
positions and slopes of the initial spectra are 
well reproduced and that the discrimination be- 
tween the five primary mass groups is sufficient. 
For scrutinizing the unfolding procedure, different 
mathematical ways of unfolding (Gold-algorithm, 
Bayes analyses, principle of maximum entropy, 
etc.) have been compared and the results are con- 
sistent ^3 . For generating the kernel functions a 
large number of EAS has been simulated [17118) 
employing CORSIKA [TS] with the hadronic in- 
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Figure 3. Results of the unfolding procedure 
using QGSJET (left) and SIBYLL (right) as 
hadronic interaction model |17) . 
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Figure 4. x 2 deviation between the forward 
folded reconstructed and measured (N e -N^)- 
data cells for the QGSJET (left) and SIBYLL 
(right) hadronic interaction models |17| . 



teraction models QGSJET (version 2001) [2U| and 
SIBYLL 2.1 HI]- 

The result of the unfolding is presented in Fig.[3| 
for each of the two interaction models. Clearly, 
there are common features but also differences in 
the energy distributions obtained with the two 
interaction models. The all-particle spectra, also 
shown in Fig. ^ coincide very nicely and in both 
cases the knee is caused by the decreasing flux 
of the light primaries, corroborating results of an 
independent analysis of Ref. [22] - Tests using dif- 
ferent data sets, different unfolding methods, etc. 
show the same behavior ^8]. As the most strik- 
ing difference, SIBYLL suggests a more promi- 



nent contribution of heavy primaries at high en- 
ergies. This difference results from the different 
N e -Nfi correlation shown in Fig.EJ i.e. SIBYLL 
predicts higher electron and lower muon numbers 
for given primaries as compared to QGSJET. 

Is there a way to judge which of the two mod- 
els describes the data better? This is done most 
easily by comparing the residuals of the unfolded 
two-dimensional N e vs AT' 1 distributions with 
the actual data used as input to the unfolding 
(Fig. |2J). The result of such an analysis is pre- 
sented in Fig.0|in terms of % 2 . The deviations re- 
veal some deficiencies of QGSJET at low electron 
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and muon numbers and they clearly demonstrate 
that SIBYLL encounters problems in describing 
the high-A^e - low-TV* 1 tail of the experimental 
data at 10 PeV and above J7|. If not being pre- 
pared to accept an additional significant contri- 
bution of superheavy primaries (A > 60) required 
in case of SIBYLL simulations to fill the gap at 
high muon numbers, the results suggest a muon 
deficit (a/o electron excess) in this model. Def- 
initely, this problem needs further attention and 
will be very important also for composition stud- 
ies at higher energies |23| . 

With this caveats kept in mind, the KASCADE 
data favor an astrophysical interpretation of the 
knee and are in agreement with a constant rigid- 
ity of the knee position for the different primaries. 
Similar results were very recently obtained for an 
analysis of two-mass groups based on combined 
EAS-TOP / MACRO measurements and 
were again confirmed for three mass groups from 
EAS-TOP electron and muon measurements |25j . 
Within the given error bars, the mean logarith- 
mic masses of both experiments agree well with 
one another. 
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Figure 5. KASCADE upper limits (95%) of 
Rayleigh amplitudes A vs primary energy (bold 
line) compared to results from other experiments 
and to expectations from galactic CR diffusion 
(thin lines) |2*B] . 



2.2. Search for Anisotropies and Point 
Sources 

Additional information about the origin of CRs 
and their propagation in the galactic environment 
can be obtained from global anisotropics in their 
arrival directions. Model calculations show that 
diffusion of CRs in the galactic magnetic field 
can result in anisotropies on a scale of 10 -4 to 
10 -2 depending on the energy of the particle, the 
strength and structure of the galactic magnetic 
field 6 , and on the source distribution. Since the 
diffusion scales again with the rigidity, a factor of 
5-10 larger anisotropies are expected for protons 
as compared to iron primaries. This rigidity de- 
pendent diffusion is one of the possible explana- 
tions of the knee. 

Because of the small anisotropy expected, a 
large data sample and careful data selection is 
necessary. About 10 s EAS events in the energy 
range from 0.7 to 6 PeV were selected and studied 
in terms of Rayleigh amplitudes A and phases $ 
of the first harmonic. Neither for the full set of 
data nor for electron-rich and -poor EAS signifi- 



cant Rayleigh amplitudes were found. The upper 
limit on the large scale anisotropy is depicted in 
Fig-El 26; an d is in line with results reported from 
other experiments. We shall come back to this re- 
sult in the next section. 

Even though the location of CR sources should 
be obscured due to the deflection of charged par- 
ticles in the magnetic field of our galaxy, there is 
interest to perform point source searches. For ex- 
ample, neutrons are not deflected and can reach 
the Earth if their energy and hence decay length 
is comparable with the distance of the source. A 
decay length of 1 kpc corresponds to a neutron 
energy of about 10 17 eV. Also, by applying appro- 
priate cuts to electron and muon numbers from 
EAS, searches for 7-ray point sources can be per- 
formed in the PeV range. 

Such a study has been performed based on 47 
Mio EAS with primary energies above ~ 300 
TeV. A certain region in the sky is then ana- 
lyzed by comparing the number of events from 
the assumed direction with an expected number 
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Figure 6. 90 % upper limit for a point source mov- 
ing through the zenith in comparison with results 
from other experiments [22]. Note the different 
definitions of the energy thresholds. 



of background events. For the latter, the so-called 
time-shuffling method has been used. As a result, 
again no significant excess is found in the region 
of the galactic plane or for selected point source 
candidates. Assuming equal power laws in the en- 
ergy spectra of background and source events, up- 
per flux limits can be calculated for given energy 
thresholds. For a steady point source that tran- 
sits the zenith, we obtain an upper flux limit of 
3 • KT 10 m^s" 1 (see Fig.EJ [21]. This is roughly 
1-2 orders of magnitude larger than the Crab flux 
extrapolated to this energy. 

Very recently, Chilingarian et al. reported the 
detection of a source of high-energy CRs in the 
Monogcm ring [2Hj- Changing slightly our cuts 
in zenith angle to widen the declination range 
thereby covering the position of the source candi- 
date, we find 742 events within an opening angle 
of 0.5° around the suggested location with an ex- 
pected number of 716 background events yielding 
an upper flux limit of 3 ■ 1CP 10 m _2 s _1 . Simi- 
lar values are found when searching for an excess 
from the direction of the pulsar PSR B0656+14 
|29| located near the centre of the Monogem SNR. 



2.3. Implications for understanding the 
CR origin 

The new high quality data presented in the pre- 
vious sections have revitalized the interest to un- 
derstand both the origin of CRs in the knee region 
and the phenomenon of the knee structure itself. 

This is because discriminating models of max- 
imum acceleration energy from galactic diffu- 
sion/drift models of the knee or from particle 
physics interpretations require detailed inspec- 
tion of knee structures seen in individual mass 
groups combined with precise measurements of 
CR anisotropics. Of particular interest are the 
energies of the spectral breaks, the power-law in- 
dices below and above the knee, and the smooth- 
ness of the turn-over regions. Even though, these 
goals are not yet achieved totally, important steps 
have been made towards it. Previous investi- 
gations were limited to inclusive CR all-particle 
spectra and to global changes of the mean log- 
arithmic mass, In A, with primary energy. Such 
measurements appeared to be too insensitive for 
a convincing discrimination of models. Further- 
more, in most cases not full distributions but only 
mean values of experimental and simulated distri- 
butions were compared to each other. A promi- 
nent example are plots of the shower maximum 
X m ax vs primary energy. Obviously, deficiencies 
of hadronic interaction models remain unrecog- 
nized in such plots, unless EAS data are below 
proton or above iron simulations. 

Very good examples demonstrating the dis- 
crimination power of the new data presented 
here and showing the amount of information con- 
tained in it is given by recent studies of Wick et 
al. [3U| or Dar [3J. Based on the earlier sug- 
gested connection between Gamma-Ray Bursts 
(GRBs) and ultrahigh-energy CRs |32I33 | they 
propose a model for the origin of CRs from 
~ f0 14 eV/nucleon up to the highest energies 
(<; f0 20 eV). In that model, GRBs are assumed 
to inject CR protons and ions into the interstellar 
medium of star-forming galaxies - including the 
Milky Way - with a power-law spectrum extend- 
ing to a maximum energy ~ 10 20 eV. High-energy 
CRs injected in the Milky Way diffuse and escape 
from our Galaxy. Ultra high-energy CRs with en- 
ergies 10 17 to 10 18 eV that have Larmor radii 
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Figure 7. Fit of the GRB-model of Wick et al. [H| 
to the preliminary KASCADE data presented at 
ICRC 2001 |34I35| . In the model, a GRB occurred 
210000 years ago at a distance of 500 pc and in- 
jected 10 52 ergs into CRs. The CRs isotropically 
diffuse with an energy-dependent mean-free-path 
in a MHD turbulence field. 



comparable to the size scale of the galactic halo 
escape directly from the Milky Way and prop- 
agate almost rectilinearly through extragalactic 
space. By the same token, UHECRs produced 
from other galaxies can enter the Milky Way to 
be detected. UHECRs formed in GRBs through- 
out the universe then travel over cosmological 
distances and have their spectrum modified by 
energy losses, so an observer in the Milky Way 
will measure a superposition of UHECRs from 
extragalactic GRBs and HECRs produced in our 
Galaxy. 

Thereby, the CR spectrum near the knee 
is understood by CRs trapped in the Galac- 
tic halo that were accelerated and injected by 
an earlier Galactic GRB. Assuming magneto- 
hydrodynamical turbulence superposed to the 
galactic magnetic field, a fit of the model of 
Ref. [301 to the preliminary KASCADE data, 
shown in Fig. suggests a 500 pc distant GRB 
that released 10 52 ergs in CRs if the GRB took 
place about 210000 yrs ago. Keeping in mind 
the still large uncertainties of the data and some 



freedom of parameters in the model, there is re- 
markable accordance observed. In this model, the 
rigidity dependence of knee position is caused by 
galactic modulation effects. The GRB-model of 
Dar 31 , on the other hand, predicts the knee to 
be proportional to the mass A of CR primaries. 
This is because of the relativistic beaming effect 
(cannon balls) in the SN-jets of that model. It 
also fits the data within the present uncertainties 
with differences showing up mostly for protons 
and Helium (Z/A — 1,0.5, respectively). 

Since in these models very few or just a sin- 
gle galactic GRB is responsible for most of the 
CRs in the knee region, anisotropics in the ar- 
rival directions are expected on different levels, 
depending on the distance and age of the GRB. 
For example, the authors of Ref. [HD] state that if 
an anisotropy below ~ 0.2 % is confirmed, then a 
number of implications follow. Either we are lo- 
cated near a rather recent GRB, which could be 
unlikely, or the CR energy release from GRBs is 
larger than the one given above [3U]. Thus, im- 
proving the anisotropy limits of the previous sec- 
tion would help to further pin down this model. 
Similarly, recent calculations of standard SN ac- 
celeration models show that global anisotropies 
well in excess even of 0.2 % are expected for more 
realistic source distributions in the Galaxy, de- 
pendent on the structure of the magnetic fields 
|3fi| . The upper limits provided by the present 
data rule out already a large parameter space in 
these calculations. 

However, before starting to over-interpret the 
data, particularly the energy spectra of mass 
groups, we should emphasize again their sensitiv- 
ity to the applied interaction models. More work 
is still needed to improve the models and to ar- 
rive at smaller systematic uncertainties. On the 
other hand, the very valuable study of such mod- 
els of CR origin demonstrates the informational 
content reached by present data. 

3. Summary and Outlook 

KASCADE has provided a wealth of new high 
quality EAS data in the knee region giving im- 
portant insight into the origin of the knee and 
of CRs in general. Conclusive evidence has been 
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reached on the knee being caused by light pri- 
maries (p + He) mostly. Furthermore, the data 
are in agreement with a rigidity scaling of the 
knee position giving support to an astrophysical 
origin by either maximum acceleration or diffu- 
sion/drift models of propagation. Astrophysical 
parameters start to be constrained by the prelim- 
inary KASCADE data, as was demonstrated at 
the example of the GRB models of Refs. |30I31| . 
Further important constrains result from mea- 
surements of large scale anisotropies of CRs and 
by new limits on point sources. Claims that were 
made in Rcf. [2E] about an excess of CRs from 
the Monogem SNR cannot be confirmed. 

A particle physics interpretation of the knee 
appears to be excluded with a high level of con- 
fidence. For example, interactions of CRs with 
background particles or photo-disintegration in 
the Galaxy would produce an abundance of sec- 
ondary protons and result in a light mass compo- 
sition above the knee energy, a result which is in 
contradiction to the present data. Furthermore, 
j7 e 's with a mass of ~ 0.5 eV/c 2 as are needed 
in the model of Ref. [H] appear to be excluded 
by recent WMAP and 2dF data for neutrinos in 
case of degenerated masses, with the latter being 
suggested by recent oscillation and possibly by a 
neutrinoless double-beta experiment (for a recent 
discussion of this topic the reader is referred to 

Refs. E3,|SH1)- 

Presently, more data and more observables are 
being analyzed within KASCADE, particularly 
in terms of composition analyses employing re- 
constructions of the muon production height |39| . 
Together with measurements of energetic hadrons 
in the central calorimeter, the unfolding tech- 
nique of electron and muon numbers in EAS has 
become a powerful tool to reconstruct the prop- 
erties of primary particles in EAS and it also 
provides important clues on how to improve the 
hadronic interaction models employed in COR- 
SIKA air shower simulations. 

KASCADE-Grande has just started its rou- 
tinely data taking and will extend the measure- 
ments up to 10 18 eV, thereby allowing to verify 
the existence of the putative Iron knee marking 
the so-called second knee in the all-particle CR 
spectrum (40141) . This, together with improved 



statistics for anisotropy measurements will allow 
to confront astro- and particle-physics motivated 
models of the knee in much more detail to the ex- 
perimental data as has been possible up to now. 

Furthermore, the use of radio antennas com- 
plementing the experimental KASCADE-Grande 
set-up may open a new window to future EAS ob- 
servations on large scales [4214 lj . Interesting first 
observations have been made and are presently 
being analyzed in detail. 
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